Vortex dynamics, pinning, and critical currents in a Ginzburg-Landau type-II 

superconductor 
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The dynamics of vortices in a type-II superconductor with 
defects are studied by solving the time-dependent Ginzburg- 
Landau equations in two and three dimensions. We show that 
vortex flux tubes are trapped by volume defects up to a criti- 
cal current density where they begin to jump between pinning 
sites along static flow channels. We study the dependence of 
the critical current on the pinning distribution and find for 
random distributions a maximum critical current equal to a 
few percent of the depairing current at a pinning density three 
times larger than the vortex line density. Whereas for a regu- 
lar triangular pinning array, the critical current is significantly 
larger when the pinning density matches the vortex line den- 
sity. 
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In a type-II superconductor, dissipation is associated 
with the motion of the vortex lattice [0J|]. This dissi- 
pation is reduced by the presence of defects, which pin 
the vortex lattice up to a critical current density where 
depinning occurs. In many applications such as super- 
conducting magnets, one is interested in optimizing the 
vortex pinning to achieve the maximum critical current. 
However, the details of the depinning transition are com- 
plex involving the non-equilibrium dynamics of an elastic 
lattice through a disordered medium. Theoretical stud- 
ies based on molecular dynamics simulations suggest the 
existence of various dynamical phases of vortex motion 
including plastic flow, uncoupled static channels and cou- 
pled channels ||,^. It is also possible to simulate vor- 
tex dynamics by solving the time-dependent Ginzburg- 
Landau equations , where the vortex- vortex interac- 
tion is completely characterised by the Ginzburg-Landau 
parameter, k. However, three dimensional Ginzburg- 
Landau vortex dynamics simulations are computation- 
ally intensive, in part because the standard explicit inte- 
gration methods require very small time-steps. In con- 
trast, semi-implicit methods are second order accurate in 
time allowing large time-steps. Although semi-implicit 
methods are widely used to simulate three dimensional 
vortex dynamics in dilute Bose-Einstein condensates ^ , 
the Ginzburg-Landau equations, involving coupled time- 
dependent vector fields, are more complex. 

In this paper we develop a semi-implicit method to 
solve the time-dependent Ginzburg-Landau equations in 
three dimensions For intermediate values of k, the 
semi-implicit method is two orders of magnitude faster 
than explicit methods, making it feasible to study dy- 



namical vortex phases, depinning, and the dependence 
of the critical current on the density and distribution of 
pinning sites. Although pinning may arise due to mag- 
netic defects, dislocations, grain boundaries, and corre- 
lated disorder such as twin planes in high-Tc supercon- 
ductors, we restrict the current study to volume defects 
which exclude the supercurrent. 

The time-dependent Ginzburg-Landau equations can 
be written as 



dtA^[yS-A)\il;[^ 



X V X A 



(1) 
(2) 



where tp is the order parameter, A is the vector potential, 
S is the phase of ij)^ and n = A/^, where A and ^ are the 
penetration depth and coherence length, respectively. In 
equations (1) and (2), distance is measured in terms of 
^, time in terms of relaxation time, r ~ /D, where D 
is the diffusion constant, and the magnetic field in terms 
of the upper critical field, Hc2- In addition, the Meiss- 
ner state critical field is given by He = 1/\/2k, and the 
depairing current density by jo — 2/3\/3 — 0.385 ||]. 
The equations are discretized using a grid of 51x51x51 
points with a grid spacing h = 0.4. The gauge invari- 
ance of the discretized equations is preserved by introduc- 
ing hnk variables of the form Uf^f. = exp(— iAfjj,/i) llo| . 
The discretized equations are solved using a semi-implicit 
Crank-Nicholson method with a time step, 6t — 0.5 A 
current flow along x is induced by imposing a magnetic 
field difference, ^Bz, between the upper {y = 10) and 
lower {y = —10) boundaries. The supercurrent across 
the boundary is set to zero. We impose periodic bound- 
ary conditions at a; = ±10 and z = ±10. The average 
current density is given by, j = hi^^Bz/d, where d is the 
width of the superconductor. A pinning array is produced 
by adding a potential term to equation (1) consisting of 
a random distribution of cubic potential steps with side 
length a = 1.2 and height Vq = 5.0. In agreement with 
other studies ||ll[, we find that the pinning strength in- 
creases with a for a < ^, and saturates for a > ^. A more 
sophisticated pinning model would be needed to account 
for the large pinning forces observed for small defects |l2| . 

Fig. 1 shows a sequence of images illustrating the mo- 
tion of the vortex lattice through the pinning array. In 
frame 1, six flux tubes are visible. By comparing frames 
1, 2 and 3, one sees that the central flirx tubes are moving 
whereas the two pairs on either side are pinned. How- 
ever, between frames 4 and 5 the flux tubes on the left 
and right jump to the next pinning site. This differential 
motion between neighbouring planes in the vortex lattice 
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plays an important role in the voltage-current character- 
istic (see below). After frame 6, a similar but not identical 
sequence recurs. For the simulations presented in Fig. 1, 
the bending of the vortex lines is increased by the choice 
of a larger value of k and strong pinning. However, no 
entangling of vortex lines is observed. For smaller k, the 
vortex lines become more rigid, and the behaviour of the 
three dimensional system and a two dimensional cross 
section are qualitatively very similar. For high- Tc su- 
perconductors, a comparison between the two and three 
dimensional dynamics should consider possible effects of 
the layered structure 
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FIG. 1. A sequence of three dimensional images showing 
the motion of a /t = 5 vortex lattice through a random pinning 
array. The axes are shown inset in frame 6. The current 
flows along a;, the external magnetic field is along z, and the 
vortices move in the —y direction. Each frame shows a region 
with dimensions 9 x 7 x 20 coherence lengths containing 12 
pinning sites (shown in black, not to scale). The external 
magnetic field and current are Bext = 0.4 and j = 5 x 10~^, 
respectively. The grey flux tubes corresponds to surfaces of 
constant supercurrent density, |?/>|^ = 0.05. The time interval 
between successive frames is 100. 

We use two dimensional simulations to study the effect 
of pinning on the voltage-current characteristic qxV — 1 
curve of a superconductor with k = 3, where three di- 
mensional effects are suppressed. In addition, we reduce 
the size and strength of the pinning sites to a = 0.8 and 
Vo — 2.0, respectively. In Fig. 2 we present contour plots 
illustrating the vortex lattice in two dimensions. Fig. 2(a) 
shows the instantaneous vortex distribution for a perfect 



superconductor (no pinning). The vortex density is pro- 
portional to the magnetic field which decreases linearly 
from the bottom to the top. The vortices move upwards 
with a speed v — E/B, where B is the local magnetic 
field and the electric field, i?, is constant throughout the 
sample. Consequently, the vortex flow obeys a Bernoulli- 
like equation where the flow is faster in regions of lower 
density (lower magnetic field) and the dissipation can be 
thought of as a relaxation of the magnetic flux lattice. 



a) 

















^^^^ 






















®« 












i 






i 




1 


Is 






4 












FIG. 2. Contour plots of the supercurrent density for a sec- 
tion of a superconductor with dimensions 60 x 20 coherence 
lengths, and k = 3, subject to an external magnetic field in 
the z direction, Bcxt ~ 0.4. The current flows in the x direc- 
tion and the vortices move in the y direction, (a) For a current 
j = 0.06 and no pinning, the vortices form a triangular lat- 
tice with lattice spacing proportional to the local magnetic 
field, (b) The addition of pinning (density 0.056 ^~^) creates 
a vortex glass, which at low currents, j = 0.02, is pinned, 
(c) At intermediate currents, j = 0.025, vortex motion begins 
along channels, indicated by the grey scale image of the local 
electric field, superimposed on a time-averaged contour plot 
of the supercurrent density, (d) At larger currents, j — 0.054, 
all the vortices are moving and the electric field is non-zero ev- 
erywhere, however, the channels, where vortex motion mainly 
occurs, are still visible. 

Adding defects transforms the triangular lattice into an 
irregular vortex glass. Fig. 2(b). For low driving fields, 
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the vortex glass is frozen. As the current is increased, in- 
dividual vortices begin to jump between pinning sites. As 
in the three dimensional simulations, Fig. 1, this motion 
begins along channels. The existence of static channels 
confirms the results of molecular dynamics simulations 
j^,^ . In the Ginzburg-Landau model channels can merge 
or divide at intermediate drive currents, as shown in Fig. 
2(c). At larger currents, all the vortices are moving but 
the channels are still evident. Fig. 2(d). 

The on-set of vortex motion coincides with the on-set of 
dissipation or breakdown of superconductivity. In Fig. 3 
we plot the V — I curve for a two-dimensional thin film 
for different defect densities. The voltage is measured 
by decreasing the current at a very slow rate of —1.2 x 
10^^ in 2.5 X 10^ steps, and the V — I curves is obtained 
from a 200 point moving average. As our sample size is 
relatively small, surface effects tends to dominate. The 
critical current due to the Bean-Livingston barrier for 
vortices entering and leaving the calculation region Q is 
the same order of magnitude as the pinning effect. In 
order to study pinning only we remove the surface effects 
by adding a boundary layer of width 9 ^ on either side of a 
calculation region with width 30 Within the boundary 
layer, a linear ramp potential reduces the supercurrent 
density gradually to zero. The current density and the 
voltage are measured within the calculation region {\y\ < 
15) only. 
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FIG. 3. The V — I curves for a two dimensional section of 
superconductor with pinning densities (from the right) 0.14, 
0.28, 0.39, and 0.56 (at these relatively high densities, 
the critical current decreases with increasing pinning density) . 
The thin black line corresponds to the V — I curve without pin- 
ning, and the dotted line shows the normal resistance, E = j. 
Note that at large currents the slope of the V — I curves is 
similar to the normal resistance curve. 

The shape of the V—I curve is dependent on the details 



of the vortex dynamics. The characteristic 'curved foot' 
can be explained by the combination of an increase in the 
number of vortex flow channels and increased flow along 
each active channel, as illustrated in Fig. 2(c) and (d). 
The V — I curve becomes linear when all the vortices 
start to move. The ratio between the V — I curves and 
the normal resistance (the dotted line in Fig. 3) gives the 
dimensionless resistivity, which measures the fraction of 
current carried by normal electrons. 

As the current is decreased the voltage becomes zero, 
i.e., all the vortices become pinned, at some finite cur- 
rent which we define as the critical current density, jc- 
In the absence of finite temperature induced fluctuations 
or vortex creep, the value of jc is well defined. However 
the critical current is sensitive to the exact distribution 
of pinning sites, therefore we average over six random 
distributions with the same density. Fig. 4 shows a plot 
of the average value of jc against pinning density. The 
maximum critical current density is about 2 % of the de- 
pairing current, jo- For comparison, the optimum critical 
current density of Nb-Ti alloy is ~ 3 % of jo ■ The max- 
imum value of jc occurs at pinning density about three 
times larger than the vortex line density (indicated by 
the dotted line in Fig. 4). 
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FIG. 4. The critical current density as a function of the 
defect density (in units of for both random distributions 
(v) and regular triangular arrays (•). The data points are 
determined from an average of six random distributions. The 
error bars (shown for the high density distribution only) in- 
dicate the standard deviation. An example illustrating the 
effect of the distribution on the V — I curves is shown inset. 
The bold curve is a fit using the function Ax exp{ — Bx), where 
A and B are fit parameters. The critical current density for a 
regular triangular array is a maximum when the pinning den- 
sity is equal to vortex line density (indicated by the dotted 
line) . 

The dependence of the critical current on the defect 
density fits reasonably well to a function of the form 
Ax eyip{—Bx). The linear increase at low pinning density 
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follows from the linear dependence of the critical current 
on the pinning force. The exponential decrease at large 
pinning densities is due to the competing effect of super- 
current depletion by defects. The shape of the curve and 
the relatively high optimum pinning density also agree 
qualitatively with experimental results on silver doped 
high-Tc superconductors JTit . 

For certain random distributions one finds persistence 
static channels which can dramatically reduce the crit- 
ical current. This is illustrated in Fig. 4(inset), where 
the curve with lower dissipation at large currents has a 
significantly lower critical current. 

One approach to increase the critical current is to intro- 
duce a regular pinning array by nanostructuring [T5|-[l^] . 
In Fig. 4 we show that a regular triangular array in- 
creases the critical currents by more than a factor of two, 
however, the optimum pinning density is sharply peaked 
around the vortex line density. Consequently, the en- 
hancement is only obtained within a narrow range of the 
external magnetic field. This agrees with experimental 
studies where a sharp enhancement peak is obtained at 
matching magnetic field values [ pT| . There are two addi- 
tional critical current peaks, one at one third the vortex 
line density where every third vortex is trapped, and one 
at half the vortex line density, which is weaker because 
the matching only occurs on alternate planes. For small 
pinning sites (a = 0.8 compared to the vortex cores size 
of 2) the maximum critical current is about 5 % of the 
depairing current, jo- For a = 2 we obtain jc = 0.074j£), 
which suggests that other pinning mechanisms may be 
needed to obtain jc ~ jo ■ 

In summary, we have studied vortex dynamics and pin- 
ning in a three dimensional superconductor by solving 
the time-dependent Ginzburg-Landau equations. We find 
that above a critical current density vortex flux tubes 
jump between pinning sites following specific channels. 
The main features of the dynamics are reproduced by two 
dimensional simulations. We study the effect of pinning 
on the voltage-current characteristic of the superconduc- 
tor, and show that the breakdown of superconductivity 
is associated with the appearance of channelled vortex 
flow. The characteristic curved foot in the V — I curve 
arises due to the combination of the formation of more 
channels and faster vortex flow along each channel. For 
a random pinning array we find a maximum critical cur- 
rent equal to 2 % of the depairing current occurring at a 
pinning density of about three times the vortex line den- 
sity. Finally, we study the critical currents produced by 
vortex matching pinning arrays. The results suggest that 
time-dependent Ginzburg-Landau simulations are ideally 
suited to provide quantitative predictions of critical cur- 
rents in type-II superconductors. 
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